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Abstract 


Singularities inside the static spherically symmetric black holes in the SU(2) 
Einstein-Yang-Mills and Einstein-Yang-Mills-dilaton theories are investi¬ 
gated. Analytical formulas are presented describing oscillatory and power 
law metric behavior near spacelike singularities in generic solutions. 
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Black holes solutions in the theories including non-Abelian massless vector fields [jl| 
are known to exhibit unusual features, which are not shared by the “classical” vacuum or 
electrovacuum black holes. Besides violation of the naive no-hair and uniqueness theorems, 
they demonstrate essentially new type of internal structure 0. It was shown that the 
space-time inside a generic static black hole in the pure Einstein-Yang-Mills (EYM) theory 
has no Cauchy horizons while the metric oscillates with an infinitely growing amplitude as the 
singularity is being approached. Although some particular solutions still may possess either 
Schwarzschild (‘S’) or Reissner-Nordstrom (‘RN’) type singularities |^, these configurations 
constitute only a zero measure set in the solution space. Inside a generic EYM black hole 
the mass function exhibits successive exponential jumps followed by abrupt falls to almost 
zero values. The amplitude of jumps is exponentially growing up towards the singularity, 
while the duration of cycles tends to zero. The behavior of the system near the singularity 
is described well by a two-dimensional dynamical system 0, which clearly demonstrates an 
infinitely oscillating regime. These results were confirmed in the existence proof for the 
RN-type solution was announced in 0. 

The purpose of the present paper is to discuss the singularity structure in generic static 
spherical black holes in the theories EYM and EYM with dilaton (EYMD) (see 0 and 
references therein). We assume the ‘stringy’ value of the dilaton coupling constant and 
without loss of generality choose the units and the length scale in which both the Planck 
mass and the gauge coupling constant are equal to unity: 

5 = / {-R + 2(V0)^ - e-^^F^} , (1) 

where F is the SU(2) field corresponding to the connection 

A^Tadx^^ = {W{r) - l){T^de - Tg sin Odcj )), 

(T^ and Tg being spherical projections of the SU(2) generators). 

It is convenient to parametrize the space-time interval as 

= — ^dt^ — T^r"^ ~ r‘^{d9‘^ + sin^ 6d(p ^), (2) 

A 


where A = — 2rm{r) is negative inside the black hole of the generic type. 

The equations of motion for W, A, </> can be decoupled from the equation for a, so that 
in the general case of the EYMD system the equations take the form 

AU' - 2AU(j)' = WV/r - TW' , (3) 

(A/r)' + A0'2 = F - 2AU‘^e-‘^^ , (4) 

(A0')' + Ar0'3 = F- 2A(0V + - 1, (5) 


where 


W' 1/2-20 

t/=—, -R=l-R 


lE^ - 1. 
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( 6 ) 


The remaining equation for a reads 

(lncT)' = r(0'2 + 2f/V2^) . 

The black hole solutions are generated by the initial data at the events horizon rh- 
hh/i = W{rh) and (ph = subject to the condition > 0. They can be labeled by 

the ADM mass M and the dilaton charge D = — hm(r^0') for r —>■ oo. The units for these 
quantities may conveniently be hxed requiring (p{oo) = 0. Then both Wh and cph are entitled 
to take the discrete values to £t the asymptotic flatness conditions. For all physical EYMD 
black holes M > D. 

Let us start with examining the singularity structure of the exponentially oscillatory 
interior space-time for generic black holes in the pure EYM theory. The corresponding 
truncated system reads 


AU'+ 1- 




W' = 


wv 


(v)' 

d 


( — ) + 2AU^ = 1 - 




dr^ 


\na = U^ 


( 7 ) 

( 8 ) 

(9) 


It is generically observed during numerical integration from the events horizon towards 
the origin (apart from some discrete horizon data) that the oscillations of A take place so 
that, when oscillation progress, the right side of the Eq. ( 0 ) becomes small with respect to 
the terms at the left side. Neglecting it, one obtains the following approximate hrst integral 
of the system: 


Z = AUa/r = const, 


( 10 ) 


which relates oscillations of the mass function to variation of a. Numerical experiments also 
show that while the YM function W remains almost constant up to r = 0, its derivative 
is still non-zero and is rapidly changing on some very small intervals of r. The function U 
exhibits a step-like behavior being constant with high accuracy during almost all the chosen 
oscillation cycle (Fig. 1) and then jumping to a greater absolute value corresponding to the 
next cycle. It is clear from (^) that a is exponentially falling down with decreasing r while 
U ~ const, whereas in the tiny intervals of D-jumps a remains almost unchanged. So a 
tends to zero through an inhnite sequence of exponential falls with increasing powers in 


the exponentially decreasing intervals. Combining this with (p!0[) and the above mentioned 
properties of U one can deduce rather detailed description of the metric behavior. 

Let us denote by r*, the value of radial coordinate where A has k-th local maximum. 
Soon after passing this point, the function U stabilizes at some value Uk approximately 
equal to the doubled value at the point of local maximum (similarly, U increases by about a 
factor of two when approaching the local maximum, whereas A is almost stationary). Then, 
according to (I), a is equal to 

a{r) = a{rk) exp \u^{r^ - rl) 


From ([10|) one hnds that, while Uk ~ const, 
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( 11 ) 


A/ A 

A(r) = -r exp 


rk 


UM-r^) 


This function falls down with decreasing r until it reaches a local minimum at 


1 \/|^(^fc)l 

* “ V2 \U,\ ~ nV(rk)\ 


( 12 ) 


In what follows, in view of the observed fact that in the course of oscillations of A the YM 
function W changes insignihcantly, we will put V = const. 

Therefore, the mass function is inflating exponentially while r decreases from to Rk- 
After passing R^, an exponential in (pd]) becomes of the order of unity, hence A starts to 
grow linearly, and the mass function m{r) stabilizes at the value = m{Rk). Such a 
behavior holds until the point of local maximum of A/r^ is reached; this takes place when 
A —V"^ at the point 


* 

r 

k ~ 


1/2 

lA/ exp 


-{UkTkf 


(13) 


After this a rapid fall of |A| is observed causing a violent rise of |17|. Then the term 2AU^ 
in the Eq. (|]) becomes negligible and consequently at this stage 

UA ^ -V^Uk , (14) 


while r practically stops. This implies that very soon A reaches the next local maximum at 
the point r^+i r^, while m(r) rapidly falls down to At the point of local maximum 

of A one has in the Eq. (|[) |A| -C then in view of the smallness of r we hnd 


\U{rk)\ 


1^1 

^2\A{rk)\rk 


(15) 


To obtain the estimates by the order of magnitude we will neglect all numerical co¬ 
efficients elsewhere except for the power indices in exponentials, in particular putting 
U{rk) = Uk, and omitting also (quasi-const ant) factors V. With this accuracy one obtains 
from ([TI|) -(p!^: 



Thus, introducing a variable Xk = {rk/RkY 1), we can derive the following recurrent 
formula 


Xk+l = , 
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which shows that Xk is an exponentially diverging seqnence. In terms of Xk one has 


Tk+l 

Tk 


Xke 




1 


this relation can also be nnderstood as a ratio of the neighboring oscillation periods since 
^k 3> Tk+i- Valnes of the fnnction |A| at the points rapidly tend to zero: 

so we deal with an inhnite seqnence of “almost” Canchy horizons as r —> 0. At the same 
time the valnes of |A| at the points Rk grow rapidly: 


|A(i?fc)| = Xfc 


— t.“ 3/2 a;fc/2 


and the valnes of the mass fnnction grow correspondingly as 

= ^-lpa:fc/2 

Mk-1 " 

As was shown in Q, in snch a regime the system (0)-(||) can be reduced to a two- 
dimensional autonomous dynamical system, one of the hxed points of which is an unstable 
focus with trajectories inhnitely spiraling outward. 

While r decreases form Vk to Rk, the function a rapidly falls down to the value ak = 
a{Rk), which then remains unchanged until the point Vk+i- As the singularity is approached, 
the sequence ak decreases according to 

Now consider the singularity structure for generic black holes in the EYMD theory. In 
this case, starting numerical integration at the horizon, one does not observe huge oscillations 
of the metric in the interior region. A generic solution does not exhibit Cauchy horizons 
either, so A remains negative for all 0 < r < r/^. For sufficiently small r the right sides of 
the Eqs. (^)-(^ become small in comparison with the left hand side terms and one gets the 
following truncated system 


(lnf/)'-20' = O, 

[ln(A/r)]' = [ln(A0')]' = • (16) 

Its integration gives the following hve-parameter (i.e. generic) family of solutions 

W = Wo + ^ 

A = , (17) 

0 = c -f In , 

with constant Wq, b, c, /i, A. The validity of the truncated equations ([T6|) now can be checked 
by substituting the asymptotic solution (pTl) into the full system (i-(l)- For consistency it 
is sufficient that the following inequalities hold: 
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\/2 - 1 < A < 1, 


which is in agreement with the numerical data. 

From (^) it follows that the mass function diverges as r —>■ 0 according to the power 
law: 

m(r) = . 

The corresponding a tends to zero as 

a(r) = , 

where ai = const. Thus, the dilaton completely suppresses the oscillations of the metric 
inside the black hole. 

Let us discuss now the geometry of space-time near the singularity. The regimes described 
above correspond to an expectation that singularity inside black holes of the general type 
should be spacelike. The metric in the interior region of a black hole corresponds to the 
anisotropic Kantowski-Sachs cosmological solution. One can show that the corresponding 
shear parameter a dehning the shear tensor aij = (2, —1, —l)a'/3, inhnitely grows as the 
singularity is approached. For an oscillating EYM solution the values of a at the points 
are of the order of 


-max 


rsj 


M, 


1/2 


p3/2 


while the values at the points are 


Mk- 

Rk- 


Both sequences are inhnitely growing. For the EYMD solutions the divergence of a is 
power-like. 

Thus, spherical non-Abelian black holes of the generic type are conformable with the 
strong cosmic censorship principle (singularity is spacelike). In the case of the pure EYM 
theory, there is, however, an inhnite sequence of “almost” Cauchy horizons in the vicinity of 
which the mass function starts to grow exponentially stabilizing on the value corresponding 
to the next oscillation cycle; the sequence of these values diverges exponentially as the 
singularity is approached. In the theory including a dilaton the mass function monotonously 
(according to the power law) tends to inhnity. In the cosmological interpretation such a 
behavior corresponds to an inhnite growth of anisotropy. 

D.V.G. thanks the Theory Division, CERN for hospitality in December 1996 while the 
work was in progress. The research was supported in part by the RFBR grants 96-02-18899, 
18126. 
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FIGURES 



FIG. 1. The first oscillation cycle for an EYM solution. The functions A for EYM (solid line) 
and EYMD (dashed) are shown in lower half-plane. In upper half-plane — mass functions m(r) 
(analogously) and the function U for EYM (dotted). All functions are power rescaled with the 
power index 1/10. Here = 4; Wh = —0.283993 for EYM, Wh = —0.298357, (t>h = 0.05623 for 
EYMD (asymptotically flat solutions with one node of IF.) 







